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OO ■ Abstract 

We consider the three-dimensional randomly diluted Ising model and study the critical behavior 
q^ , of the static and dynamic spin-spin correlation functions (static and dynamic structure factors) at 
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the paramagnetic-ferromagnetic transition in the high-temperature phase. We consider a purely 
relaxational dynamics without conservation laws, the so-called model A. We present Monte Carlo 
simulations and perturbative field-theoretical calculations. While the critical behavior of the static 



'^ \ structure factor is quite similar to that occurring in pure Ising systems, the dynamic structure factor 

shows a substantially different critical behavior. In particular, the dynamic correlation function 

shows a large-time decay rate which is momentum independent. This effect is not related to the 
(N 
^ \ presence of the Griffiths tail, which is expected to be irrelevant in the critical limit, but rather to 



the breaking of translational invariance, which occurs for any sample and which, at the critical 
point, is not recovered even after the disorder average. 



P: ' PACS numbers: 64.60.F-, 75.10.Nr, 75.40. Gb, 75.40.Mg 



% 



I. INTRODUCTION AND SUMMARY 

The effect of disorder on magnetic systems remains, after decades of investigation, a not 
fully understood subject. It is then natural to investigate relatively simple models, to try 
to understand the common features of disordered systems. In this regard, randomly diluted 
spin systems are quite interesting. First, they represent simple models which describe the 
universal properties of the paramagnetic-ferromagnetic transition in uniaxial antiferromag- 
nets with impurities^ and, in general, the order-disorder transition in Ising systems in the 
presence of uncorrelated local dilution. Second, they give the opportunity for investigat- 
ing general problems concerning the effects of disorder on the critical behavior. Indeed, 
several important results, theoretical developments, and approximation schemes found for 
these models have been later generalized to more complex systems like spin glasses, quantum 
disordered spin models, etc.... 

In this paper we consider three-dimensional randomly diluted Ising (RDI) systems. Their 
critical behavior has been extensively studied. ^'^'^ There is now ample evidence that the mag- 
netic transition in these systems, if it is continuous, belongs to a unique universality class, 
and several universal properties are now known quite accurately. Beside the static critical 
behavior, we also investigate the dynamic critical behavior, considering a purely relaxation 
dynamics without conservation laws, the so-called model A,^ which is appropriate for uni- 
axial antiferromagnets. We focus on the dynamic (time-dependent) spin-spin correlation 
function 



G{x2 -Xi,t2- h) = {a{xi,ti) a{x2, ta)), (1) 

where a{x, t) is an Ising variable, the overline indicates the quenched average over the 
disorder probability distribution, and (■ ■ ■) indicates the thermal average. From the function 
G{x,t), one obtains the static (equal-time) structure factor G{k) and the dynamic structure 
factor G{k,uj). They are physically relevant quantities, which can be measured in neutron 
or X-ray scattering experiments.^ It is therefore interesting to study the effects of disorder 
on these physical quantities, and check whether disorder gives rise to qualitative changes 
with respect to pure systems. We investigate their scaling behavior close to the magnetic 
transition for T -^ T+ in the high-temperature phase. As we shall see, while the critical 
behavior of the static structure factor is very similar to that in pure Ising systems, the 



critical behavior of the dynamic structure factor is significantly different; in particular, the 
large-momentum behavior shows some new features. 

Since the critical region in the paramagnetic phase, i.e. for T > T^, is located in the 
Griffiths phase, ^''^ it is mandatory to discuss first the relevance of the so-called Griffiths 
singularities and Griffiths tails for the universal critical behavior of the correlation func- 
tions when T -^ T^ . In fact, one of the most notable features of randomly diluted spin 
systems is the existence of the so-called Griffiths phase for T^ < T < Tp, where Tp is the 
critical temperature of the pure system. This is essentially related to the fact that, in the 
presence of disorder, the critical temperature T^ is lower than Tp, and therefore, in the tem- 
perature interval T^ < T < Tp, there is a nonvanishing probability to ffnd compact clusters 
without vacancies (Griffiths islands) that are fully magnetized. They give rise to essential 
nonanalyticities in thermodynamic quantities. ^'^ Moreover, these clusters are responsible for 
a nonexponential tail in dynamic correlation functions. ^'^'^°'^^ In the case of RDI systems 
one can show that 

G{x,t) ^ Gait) = 5exp[-C(lnt)3/2], (2) 

for any finite x and t -^ oo, which implies a diverging relaxation time. We should men- 
tion that these effects are quite difficult to detect, and there is still no consensus on their 
experimental evidence even in systems with correlated disorder, in which these effects are 
magnified (see, e.g., Ref. 7,12 and references therein). 

Griffiths essential singularities give quantitatively negligible effects on thermodynamic 
quantities and on the static critical behavior. One can argue that also the Griffiths tail (2) 
is irrelevant in the critical limit: the nonexponential tail does not contribute to the critical 
scaling function associated with G{x,t).^ This is essentially due to the fact that B and G 
that appear in Eq. (2) are expected to be smooth functions of the temperature, approaching 
finite constants as T ^ T^. Thus, in the critical limit, t -^ oo, T ^ T^ at fixed tC,^^, 
where ^ ~ (T — T^)^'^ is the diverging correlation length, the nonexponential contribution 
simply vanishes. To understand why, let us consider the simplified situation in which the 
contributions to the autocorrelation function G{x,t) due to the Griffiths islands and to the 
critical modes just sum as 

G{t) ^ Gcit) + Gcit) = ae^exp(-ctr') + Bexp{-G{lntf/^), (3) 

where we neglect all couplings between Griffiths and critical modes. Here, Gc{t) is the 



critical contribution, while Gait) is the nonexponential Griffiths tail, which dominates for 
t ^ t*, where t* is the time at which the two terms have the same magnitude. In the critical 
limit we have t* ~ ^^(In^)^^'^. Since the critical limit is taken at fixed t/,^^, the relevant 
quantity is t* /^^, which diverges as (In,^)^/^ approaching the phase transition. This means 
that, for any fixed value of t/C,^, the condition t <^ t* is always satisfied sufficiently close to 
the critical temperature Tc, i.e., the nonexponential tail is negligible. 

In order to determine the scaling behavior of the spin-spin correlation function (1), we 
perform Monte Carlo (MC) simulations of a three-dimensional RDI model on a simple 
cubic lattice and perturbative field-theoretical (FT) calculations. In the following we briefiy 
summarize our main results. 

The high-temperature critical behavior of the static structure factor G{k) is substantially 
analogous to that of G{k) in pure Ising systems. ^^ We consider the universal scaling function 
g{Q'^) = G{k)/G{0), where Q^ = k"^^"^ and ^ is the second-moment correlation length. We 
find that g{Q^) is very well approximated by the Ornstein-Zernike form (Gaussian free-field 
propagator) g{Q)oz = 1/(1 + Q"^)'- deviations are less than 1% for Q ^ 5 and increase to 
5% at Q ~ 50. At large momenta, for <5 ^ 10 say, the static structure factor follows the 
Fisher-Langer law;^^ in particular, g^Q"^) ^ 0.92/Q'^''^ with rj ^ 0.036 for Q > 30. 

At variance with the static case, the dynamic structure factor displays substantial differ- 
ences with respect to the pure case, even though, as expected, the Griffiths tail turns out to 
be irrelevant in the critical limit. We consider the universal scaling function 

^(Q^5)^limSM (4) 

T^T+ G{k, 0) 

in the critical limit t -^ oo. A; ^ 0, and T -^ T+ at fixed Q and S. Here S = t/rjnt, where 
Tint is the zero-momentum integrated autocorrelation time. Perturbative field theory shows 
that r(Q^, S) decays exponentially, as in the case of pure systems, and this is confirmed by 
the simulation results. In pure Ising systems^^ the large-S* decay rate of T{Q^, S) depends 
on Q: the large-S* behavior is very similar to that of T{Q^, S) in the noninteracting Gaussian 
model, i.e. T{Q'^,S) ~ exp[—K{Q'^)S], where k^Q"^) = (1 + Q^). This behavior drastically 
changes in the presence of random impurities; in particular, the large-S* decay rate becomes 
independent of Q. MC simulations and FT perturbative calculations show that, for generic 
values of Q, r{Q'^,S) has two different behaviors as a function of S. For small values of 
S, it decreases rapidly, with a rate that increases as Q increases, as it does in the pure 



system. ^^ For large S instead, T{Q^, S), and therefore G{k,t), decreases with a momentum- 
independent rate. For large Q and S we find 

^(Q^S)~5'^Q-'^e-"^ (5) 

where a and ( are critical exponents, and k, does not depend on Q. We present a physical 
argument which relates the different large-S* behavior compared to pure systems to the loss 
of translational invariance. Such a phenomenon is obvious for a given fixed sample, but at 
the critical point translational invariance is not even recovered after averaging over disorder, 
because of the absence of self-averaging. Note that this phenomenon is only related to 
disorder and thus, it is expected in all systems in which disorder is relevant. 

In general, the perturbative calculations predict a scaling behavior of the form 

T{Q',S)r^Q-<f^{S), (6) 

for large Q, where /^(S*) is a function of S such that 7^(5') ~ S"'e^'^^ for S -^ oo. This 
behavior implies that G{x,t) is always nonanalytic for a; = and any t. Indeed, because of 
Eq. (6), the integral 

f d^Q Q^'T{Q\ S)g{Q^) ~ f d^k k^^'Gik, t) (7) 

diverges for n > nc = (C + 2 — rf — r7)/2 (rf is the spatial dimension). Since the moments 
of G{k,t) are directly related to the derivatives of G{x,t) with respect to x computed for 
a; = 0, the n-th derivative of G{x,t) diverges for n > ric, hence G{x,t) is not analytic at 
X = 0. This implies that the scaling function -F(V^, S), defined as 

G{x, t) = ^-^+^-^F{Y\ S), Y'^ x^e, (8) 

behaves as F{Y\ S) = fo{S) + fx{S)\Y\^ + ■■■, where X = ( + 2 - d - r], for Y^ ^ 0. MC 
simulations indicate that (^ ^ 2 in three dimensions, which implies A ~ 1. This phenomenon 
does not occur in pure systems, since in this case the decay rate k depends on Q and 
guarantees the integrability of the integrand which appears in Eq. (7) for any n. Therefore, 
F{Y\ S) has an analytic expansion around Y^ = 0, i.e., F{Y\ S) = fo{S) + f2{S)Y^ + ■■-. 
This nonanalyticity should be a general property of random models in which disorder is 
relevant, and not specific of RDI systems. 



Finally, it is worth mentioning that in the case of three-dimensional randomly dilute 
multicomponent spin models, such as the XY and the Heisenberg model, the effects of 
disorder found in RDI models are expected to be suppressed in the critical limit T — > T+, 
and should only appear as peculiar scaling corrections.^^ Indeed, the asymptotic critical 
behavior of the correlation functions is expected to be the same as that in the corresponding 
pure model, because the pure fixed point is stable under random dilution (according to the 
Harris criterion, ^^ dilution is irrelevant if the specific- heat exponent a of the pure system is 
negative). 

The paper is organized as follows. In Sec. II we introduce the model that we study. In 
Sec. Ill we discuss the static structure factor in the high-temperature phase. In Sec. Ill A we 
define the quantities that are computed in the MC simulation, in Sec. Ill B we present some 
perturbative calculations, while in Sec. Ill C we discuss the MC results. In Sec. IV we discuss 
the dynamic structure factor. Again, we first define the basic quantities (Sec. IV A), then 
we present a one-loop perturbative calculation (Sec. IV B), and finally we report the MC 
results (Sec. IV C and IV D). In the appendices we report some details of the perturbative 
calculations. 

II. THE MODEL 

We consider the randomly site-diluted Ising model with Hamiltonian 

<xy> 

where the sum is extended over all pairs of nearest-neighbor sites of a simple cubic lattice, 
(Ta: = ±1 are Ising spin variables, and p^ are uncorrelated quenched random variables, 
which are equal to 1 with probability "p (the spin concentration) and with probability 
1 — p (the impurity concentration). For Ps < p < 1, where Ps is the site-percolation point 
{Ps = 0.3116081(13) on a simple cubic lattice^^), the model has a continuous transition with 
a ferromagnetic low-temperature phase. 

MC simulations of RDI systems have shown rather conclusively (see, e.g., Refs. 1,2,3, 
19,20,21) that their continuous transitions belong to a single universality class. The RDI 
universality class has been extensively studied by using FT methods and MC simulations. 
At present, the most accurate estimates of the critical exponents are^*^ z/ = 0.683(2) and 

6 



1] = 0.036(1), obtained by a finite-size analysis of MC data. Tliese estimates are in good 
agreement witli tliose obtained by using field theory. An analysis of the six-loop perturbative 
expansions in the three-dimensional massive zero-momentum scheme gives^^ u = 0.678(10) 
and 1] = 0.030(3). Note the good agreement between FT and MC results, in spite of the 
fact that the perturbative FT series for dilute systems are not Borel summable.^^'^'*'^^ Also 
the correction-to-scaling exponents have been determined quite accurately. For the leading 
exponent u, MC simulations give^^ u = 0.29(2) (older simulations gave u = 0.33(3)^° 
and uj = 0.37(6)^9), while field theory predicts uj = 0.25(10),22 uj = 0.32(6).^^ For the 
next-to-leading exponent uj2, an appropriate analysis of the FT expansions gives^° 002 = 
0.82(8), which is consistent with the MC results for improved models. ^°'^^ Beside the critical 
exponents, also the equation of state, ^^ some amplitude ratios, ^^'^^'^'^ the universal crossover 
functions between the pure and the RDI fixed point^^ and between the Gaussian and the 
RDI fixed point, ^^'^^ and the crossover exponent in the presence of a weak random magnetic 
field^^ have been computed. 

A full characterization of the types of disorder that lead to a transition in the RDI 
universality class is still lacking. For instance, RDI transitions also occur in systems that are 
not ferromagnetic: this is the case of the Edwards- Anderson model (±J Ising model), which 
is frustrated for any amount of disorder. ^^ Nonetheless, the paramagnetic-ferromagnetic 
transition line that starts at the pure Ising transition point and ends at the multicritical 
Nishimori point belongs to the RDI universality class. ^^ 

Beside the static behavior, we also consider the critical behavior of a purely relaxation 
dynamics without conservation laws, the so-called model A, as appropriate for uniaxial 
magnets.^ The critical behavior of the model-A dynamics for RDI systems has been re- 
cently studied numerically (for a critical review of the existing results, see Ref. 21). It has 
been shown that the relaxational dynamics belongs to a single dynamic universality class, ^^ 
characterized by the dynamic critical exponent z = 2.35(2). 

For Hamiltonian (9) an accurate study of the dependence of the size of the corrections to 
scaling on p is reported in Ref. 20. It turns out that the leading scaling corrections associated 
with UJ = 0.29(2) are suppressed for p = p* = 0.800(5), in agreement with the findings of 
Refs. 19,27. For this reason we have performed our simulations at p = 0.8. 



III. STATIC STRUCTURE FACTOR IN THE HIGH-TEMPERATURE PHASE 

A. Definitions 

We consider the static (equal-time) two-point correlation function G{x) = G{x,t = 0). 
In the infinite-volume limit we define the second-moment correlation length C, 



^2 ^ _ 1 dGik) 



(10) 

fc2=0 



X dk? 
where G{k) is the Fourier transform of G{x) and 

X = 5^G(x) = G(0) (11) 

X 

is the magnetic susceptibility. It is also possible to define an exponential correlation length 
^exp- Given the infinite- volume G(x), we define 

Cexp^- lim -^L. (12) 

|x|-*oo lnG(a;) 

In the critical limit ^ and ^cxp diverge. If t^ = {T — Tc)/Tc and T^ is the critical temperature, 
for \t.f\ ^0 we have in the thermodynamic limit 

'Ci'Ccxp ~ {tA'"-, (13) 

where z/ is a universal critical exponent. In the same limit, correlation functions have a 
universal behavior. For instance, the infinite-volume G{k)/x becomes a universal function 
of the scaling variable 

Q' = k^: (14) 

i.e. we can write in the scaling limit A; — > 0, t^. — ^ at fixed Q 

X-'G{k)^g{Q'), (15) 

where g{x) is universal. Moreover, the ratio C^/^cxp converges to a universal constant Sm 
defined by 

eve ^ Sm- (16) 

For a Gaussian theory the spin-spin correlation function shows the Ornstein-Zernike (OZ) 
behavior 

Gozik) = ^. (17) 



It follows X = Z:/r^ ^"^ = 1/r, and 

gozm = Y^ ■ (18) 

Moreover, Sm = 1- 

Fluctuations change this behavior. For small Q^, g{Q'^) is analytic, so that we can write 
the expansion 

(7(Q2)-i = 1 + Q^ + 5^c„Q^ (19) 

n=2 

where the coefficients c„ parametrize the deviations from the OZ behavior. For large Q^, 
the structure factor behaves as 

12\ ^1 Ml '-^2 C's 



a behavior predicted theoretically by Fisher and Langer^'* and proved in the FT framework 
in Refs. 35,36. 



B. Field-theory results 

We determined the coefficients c„ by using two different FT approaches: the y^-expansion 
approach, in which the renormalization-group parameters are computed as series in powers 
of y^, e = 4 — d, and the massive zero-momentum (MZM) approach, in which one works 
directly in three dimensions. We computed the ffist few coefficients c„ to 0(e^'^) in the 
^/e expansion, and to four loops in the MZM scheme. The corresponding expansions are 
reported in App. A. Setting e = 1 in the y/e expansions (A6), we obtain C2 = —4 x 10"'^, 
C3 = 1.0 X 10"^, C4 = —4 X 10~^, and C5 = 2 x 10^^. In the MZM approach, resumming the 
perturbative expansions (A9) as discussed in Ref. 22, we obtain 

C2 = -4(1) X 10-^ (21) 

C3 = 1.2(3) X 10-^ (22) 

C4 = -5(2) X 10"^. (23) 

These results are fully consistent with those obtained in the y/e expansion, in spite of the 
fact that in that case we have not applied any resummation and we have simply set e = 1. 
We can also compute Sm- Since the coefficients c„ are very small, we obtain 

^M~l + C2 = 0.9996(1), (24) 



where we used the estimate (21) of C2. As in the Ising case/^'^^ the coefficients c„ show the 
pattern 

|c„| < |c„„i| < ... < IC2I < 1 for n > 3. (25) 

This is consistent with the expected analyticity properties of G{k). Since the complex-plane 
singularity in G'(A;)~^ that is closest to the origin is expected to be the three-particle cut 
located at A; = ±3i/C,cxp,^^'^^ the function g{Q'^)'^ is analytic up to Q^ = —9Sm- It follows 
that c„ ~ — c„_i/(3v^Sm), at least asymptotically. 

The large-Q behavior can be investigated in the ^/e expansion. The three-loop calculation 
of the two-point function reported in App. A 1 allows us to determine the perturbative 
expansion of the coefficients Ci appearing in Eq. (20). Setting e = 1 in the expressions (A8), 
we obtain Ci ^ 0.95, C2 + Q = -0.96. 

In order to compare with the experimental and numerical data it is important to deter- 
mine g{Q'^) for all values of Q. For the pure Ising structure factor, several interpolations 
have been proposed with the correct large- and small-Q behavior. ^^'^^'^^'^°'^^'^^'^^ The most 
successful one is due to Bray,^^ which incorporates the expected singularity structure of 
g{Q'^). In this approach, one assumes 1/ g{Q'^) to be well-defined in the complex Q^ plane, 
with a cut on the negative real Q^ axis, starting at the three-particle cut Q^ = — r^ with 
r^ = 95*^- Then, one obtains the spectral representation 

.i-n nn) 



poo 

H{Q^) = / duu^-"^ 

J r 



g2 + „2' 



1 1 I Q'Hm .2g. 

g{Q^) ^ SMHi-Suy ^ ' 

where F{u) is the spectral function, which must satisfy -F(-l-oo) = 1, F{u) = for u < r, 
and F{u) > for -u > r. 

In order to obtain an approximation one must specify F{u). Bray^^ proposed to use a 
spectral function that gives exactly the Fisher-Langer asymptotic behavior, i.e. 



where 



C2 7rp C3 71 

PAu) = I -\ cos 1 — 7^ cos — , 

T^ ( \ C2.TTpC3.Tr 

"2 M = — sm — 7^ sm — , 28 

^ ^ uP 2 mV^ 2z/' ^ ^ 
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with p = (1 — a)/v. To obtain a numerical expression we fix z/ = 0.683, t] = 0.036,^*^ and 
use the estimate (24) of Sm- We must also fix C2 and C3. Bray proposes to fix C2 + C3 to 
its e-expansion value (in our case C2 + C3 = —0.96) and then to determine these constants 
by requiring Fb{u = r) = 0. These conditions give C2 = —8.04 and C3 = 7.07. As a check, 
we can compare the estimates of c„ and Ci obtained by using Bray's approximation gsiQ'^) 
with the previously quoted results. We obtain 

C, = '^-^^^^^SmH{-Sm) ^ 0.92, (29) 

71 

and C2 ~ —4 ■ 10~^, C3 ~ 9 ■ 10~^, C4 ~ —4 ■ 10^'^. These results are in very good agreement 
with those obtained before. 

C. Monte Carlo results 

In this section we study Hamiltonian (9) at p = 0.8 in the high-temperature phase, with 
the purpose of determining the infinite- volume spin-spin correlation function G{k). We per- 
form simulations on lattices of size 32 < L < 256 in the range 0.275 < f3 < 0.2856 [note 
that^^ f3c = 0.2857431(3)]. The number of samples varies with N, being of the order of 
3000, 10000, 30000, 40000 for L = 256, 128, 64, 32. For each sample, we start from a random 
configuration, run 1000 Swendsen-Wang and 1000 Metropolis iterations for thermalization, 
and then perform 2000 Swendsen-Wang sweeps. At each iteration we measure the correla- 
tion function G{x;f3,L) and the structure factor G{k; f3,L). Since rotational invariance is 
recovered in the critical limit, to speed up the Fourier transforms, we determine it as 

G{k- f3,L) = ^ 5^(e^'=^ + e^'y + e''nHO,0,0)a{x,y,z)), (30) 

x,y,z 

where the sum runs over the coordinates {x, y, z) of the lattice sites. Of course, on a finite 
lattice k can only assume the values 27rn/L, where n is an integer such that < n < L — 1. 
We also compute the second-moment correlation length ^(/3, L) defined by 

where /cmin = Stt/L, k = 2sinA;/2. For L — ^ 00, ^{P,L) converges to the infinite-volume 
definition (10) with L~^ corrections. 
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FIG. 1: (Color online) Estimates of h{k;P,L) vs C{P,L)/L at fixed Q = W (left) and Q = 10 
(right). We only report data corresponding to A; < Ajmax = tt/S. We also report an interpolating 
curve (dashed line), which is obtained by fitting all data at fixed Q reported on the right, as 
explained in the text. 

In order to determine g{Q'^) we go through several different steps. First, for each (3 and L 
we interpolate the numerical data in order to obtain G{k\ f3, L) for any k in the range [0, vr]. 
For this purpose we fit the numerical results for h{k] (3, L) = G{0; l3, L)/G{k] (3, L) to 

/i(fc; /?, L) = 1 + ^ a„P", A; = 2 sin-. (32) 

n=l 

We increase nmax until the sum of the residuals (x^) is less than half of the fitted points 
(those corresponding to 1 < n < L — 1), i.e. x^ < L/2 (note that the data are strongly 
correlated and thus it makes no sense to require x^/DOF ^ 1, where DOF = L — 1 — nmax 
is the number of degrees of freedom of the fit). In most of the cases we take nmax = 5, but 
in a few cases we had to take n^ax as large as 10. 

Then, we investigate the finite-size effects. In the critical limit we expect 

h{k- ^,L)^f(q^ ka(3, L), ^^"j . (33) 



Equivalently, one can also use 



h{k; P,L)^F(Q^ kiP, L), %^ ) , (34) 



where k = 2 sin k/2. The two scaling forms are equivalent in the scaling limit k —>■ 0, L ^ oo, 
^{P,L) — > oo at fixed Q (or Q) and ^{P,L)/L; as a consequence, the function F{x,y) is 
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the same in the two cases. Indeed, k = k + 0{k^), and thus, by keeping fixed Q or Q, one 

only changes analytic corrections decaying as L"^. In particular, whatever choice is made, 

the structure factor g{Q'^) is equal to 1/-F((5,0). Apparently, the corrections we are talking 

about here are less relevant than the nonanalytic corrections that should decay as L"'^^, 

002 = 0.82(8), and thus, a priori one would expect only small differences between the two 

approaches. Instead, as we show below, only by keeping Q fixed is one able to determine 

the structure factor in the infinite- volume limit. 

In Fig. 1 we show the numerical data for Q = 10 (left) and Q = 10 (right). On the left 

one observes very large size corrections which make impossible in practice the determination 

of the infinite-volume limit C,{f3,L)/L -^ 0. On the right instead, there are no significant 

scaling corrections and all data fall approximately on a single curve. Size corrections are 

small for ^{(3,L)/L < 0.20 and the extrapolation to C,{(3,L)/L ^ is feasible. In the two 

panels we also show the interpolation of the data at fixed Q. As expected, the data at fixed 

Q converge to this interpolation, but it is clear that no real information could have been 

obtained on the infinite-volume limit from the data in the left panel. In order to clarify why 

scaling at fixed Q is so much better than scaling at fixed Q, we consider the lattice Gaussian 

model with nearest-neighbor couplings. In this case, the spin-spin correlation function on a 

finite lattice is given by 

Gcik) = y-^, (35) 

k^ + r 

so that ^'^ = 1/r and 

h{k) = 1 + Q2. (36) 

Thus, if we take the finite-size scaling limit at fixed Q there are no finite-size corrections: 
the scaling is exact on any finite lattice. On the other hand, at fixed Q we obtain 

1 Q^iye) , 



h{k)^{l + Q' 



1- 



(37) 



12L2 1 + g2 

In this case we have l/L^ corrections, which diverge as ^/L -^ 0, exactly as we observe 
in our data. These corrections moreover increase with Q and thus make it difficult, if not 
impossible, to estimate the structure factor. 

As a consequence of the above- reported discussion we consider below the finite-size scaling 
limit at fixed Q. In Fig. 2 we show h{k; f3, L) for Q = 5, 20, 50. Since k can be at most 2, for 
each Q we can only consider values of j3 and L such that .^(/3, L) > Q/2. However, since the 
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FIG. 2: (Color online) MC results for h{k;P,L) vs S,{P,L)/L at fixed Q for three different values 
of Q: 5 (top), 20 (middle), 50 (bottom). Only data satisfying k < k^^x = 7r/3 are reported. The 
interpolation (dashed line) corresponds to a fit of the data with L > 64 as described in the text. 
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critical limit is obtained for A; —^ 0, results close to the antiferromagnetic point k = tt cannot 
have a good scaling behavior. Therefore, in the analysis we have only considered values of 
k such that k < /Cmax- If ^max varics between 7r/4 and tt/3, the final results are essentially 
independent of k^^.^. The data reported in Figs. 1 and 2 scale as predicted by Eq. (34). 
Within the precision of our results some corrections to scaling are only visible for Q = 5 and 
C,{f3,L)/L > 0.2. They however die out fast in the interesting limit C,{l3,L)/L -^ 0. Note 
also that, as Q increases, the number of available points decreases and indeed we are not 
able to go beyond Q ~ 50 with our data. 

In order to determine the infinite- volume limit F{Q, 0), we have taken all data satisfying 
C,{f3, L)/ L < 0.5 and we have fitted them to 

J max 

h{k- (3, L)\^^^Q = ao + J2 «J- exp[-jL/e(/5, L)]. (38) 

The fitting form (38) is motivated by theory, which predicts exponentially small finite-size 
corrections in the high-temperature phase. With the precision of our data it is sufficient 
to take jmax = 2 to obtain x^/DOF < 1. The coefficient ao allows us to estimate g{Q^): 
9{Q^) = l/-^('5)0) = l/cto- The results for k^^^x = 7r/4 and n/S are essentially identical 
within errors up to Q ~ 40 (for /cmax = 7r/4 we do not have enough data to determine 
reliably F{Q,0) for Q > 40). In the following we take those corresponding to k^^ax = 7r/3, 
which allow us to compute giQ"^) up to Q = 50. In order to detect scaling corrections we 
have repeated the analysis including each time only data such that L > Lmin. The results 
are essentially independent of Lmin- For instance, for Q = 5, one of the values we considered 
in Fig. 2 (in this case some scaling corrections are present for ^{j3,L)/L > 0.20), we obtain 
ao = 25.881(6), 25.881(7), 25.882(9) for L^^^ = 32,64, 128, respectively. This is due to the 
fact that Oq is determined by the results at small values of ^(/9, L)/L and in this range there 
are essentially no scaling corrections. In the following we choose conservatively Lmin = 64. 
Our final estimate of g{Q'^) is reported in Fig. 3. Deviations from the OZ behavior are 
quantitatively small and indeed at Q = 50 the relative deviation is only 0.05. It is important 
to note that the estimates oi g{Q^) at different values of Q are correlated since the estimates 
of G{k\ (3, L) for different values of k are statistically correlated. This explains the regularity 
of the results. Note also that the error changes rather abruptly in a few cases. For instance, 
this occurs between Q = 40 and Q = 4:1. This happens because at Q = 40, the estimate of oq 
is essentially determined by the result obtained for (3 = 0.2853, L = 256, which corresponds 
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FIG. 3: (Color online) Estimates of the scaling function g(Q'^){l + Q'^) for integer values of Q. 
We also report Bray's approximation, in which C2 + C3 is fixed to the \/e value (Bray-e), and the 
structure factor for the pure Ising model (Ising). The curve "interp" (solid line) corresponds to 
the interpolation gmtiQ'^) reported in Eq. (40). 

to ^{f3,L)/L = 0.158 and k = 0.989. For Q = 41 this lattice is no longer considered, since 
the corresponding k exceeds fcmax = 7i"/3 (A;max = 1)- For Q = 41, the result with the smallest 
^(/3, L)/L corresponds to C,{l3, L)/L = 0.236. The extrapolation to the infinite-volume limit 
is therefore much more imprecise. 

For Q ^ 00, g{Q^) ^ Ci/Q^-'', see Eq. (20). We fit the estimates of \ng{Q^) reported 
in Fig. 3 (they correspond to integer values of Q between 1 and 50) to a + (r/ — 2) InQ. 
If we include only data with Q > Qmm = 15 and 20, we obtain 77 = 0.032(1), 0.032(2), 
respectively. The error we quote here assumes that all data are independent, which is not 
the case. In order to determine the correct error bar, one should take into account the 
covariance among the results at different values of Q. This is not easy and therefore, in 
order to estimate the role of the statistical correlations, we use a more phenomenological 
approach, li g^stiQ^) is the estimate oi gi^Q"^) and <j{Q'^) the corresponding error, we consider 
new data gestiQ"^) — o"(Q^) with the same error and we repeat the fit. We obtain rj = 0.029 
and T] = 0.027 for Qmm = 15 and 20. Analogously, if we consider gestiQ"^) + (^{Q"^), we 
obtain rj = 0.035, 0.037. This simple analysis indicates that ±0.005 is a plausible estimate 
of the statistical error. Therefore, we quote rj = 0.032(5) as our final result. This estimate 
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is in good agreement with that reported in Ref. 20, rj = 0.036(1), obtained from a finite-size 
scahng analysis of the susceptibihty. In order to estimate Ci, we consider g{Q'^)Q'^'^ , fixing 
r]tor] = 0.036(1).20 For Q > 20 this quantity is essentially constant: g{Q^)Q^-'' = 0.921(1), 
0.920(1), 0.917(2), 0.919(3), for Q = 20,25,30,35. We thus take 

Ci = 0.919(3) [3] (39) 

as our final estimate. The error in brackets gives the variation of the estimate as rj varies 
by one error bar (±0.001). This estimate is close to the FT result Ci ~ 0.95 and in perfect 
agreement with the estimate (29) obtained by using Bray's approximation for the spectral 
function, Ci ~ 0.92. Indeed, as can be seen in Fig. 3, Bray's interpolation represents a very 
good approximation of the numerical data, deviations being quite tiny. 

In Fig. 3 we also report the structure factor in pure Ising systems (we use the phenomeno- 
logical approximation reported in Ref. 13, see their Eq. (30) with Qmax = 15 and nmax = 6). 
In the pure case, deviations from the OZ behavior are larger: the addition of impurities has 
the effect of reducing the deviations from the OZ behavior. 

Finally, we report a phenomenological interpolation which reproduces well our numerical 
data and is consistent with the large Q'^ behavior, g{Q'^) ^ 0.919(5°'°^^/(1 + Q^): 

,^2x (1 + 0.0227953g2 + 0.0000839355g^)°°°9 

^int(g') = ^ y:^2 — • (40) 

IV. DYNAMIC STRUCTURE FACTOR IN THE HIGH-TEMPERATURE PHASE 

In this section we consider the dynamic behavior of the Metropolis algorithm, which is 
a particular example of a relaxational dynamics without conservation laws, the so-called 
model A, as appropriate for magnetic systems. In Ref. 21 we computed the dynamic critical 
exponent, obtaining z = 2.35(2). Here, we focus on the dynamic structure factor. 

A. Definitions 

To investigate the dynamic behavior we consider the time-dependent two-point correlation 
function (1) and its Fourier transform G{k, t) with respect to the x variable. Then, we define 
the integrated autocorrelation time 

1 y, G(M) _ 1 v^ G(M) ,^j, 
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and the exponential autocorrelation time 

rexp(A;) = - lim l^j (42) 

l*l-»°o \nG{k,t) 

which controls the large-t behavior of G{k,t). Here t is the Metropolis time and one time 
unit corresponds to a complete lattice sweep. 

Beside Ti^t{k) and Tcxp{k) we also define autocorrelation times T[^t,x and Tcxp,^-^^ In general, 
given an autocorrelation function A{t) we define 

Tl 

""^''^^^ " ln[A(s-n/2)/A(s + n/2)]' ^^^^ 

for any integer s and any fixed even n. By linear interpolation these functions can be 
extended to any real s. Then, we define Ti^t,x and r^^p^x as the solutions of the consistency 
equations 

'exp,x 'cS\-'^'cxp,x)j V / 

^int.x = IixTint,x)- (46) 

These definitions have been discussed in Ref. 21. There, it was shown that they provide 
effective autocorrelation times with the correct critical behavior. For a; ^ oo, Tcxp,x and 
Tint,a; couvcrge to T^^p and Tint, respectively. 

As discussed in the introduction, for Tc < T < Tp the correlation function G{x, t) does not 
decay exponentially for any finite value of x, but presents a slowly decaying tail, cf. Eq. (2). 
Therefore, Tcy^p{k) diverges for all T^ <T < Tp. As discussed in Ref. 21, this is not the case 
for the effective exponential autocorrelation time Texp,x, which is finite for any finite x. Note 
that correlation functions decaying as in Eq. (2) have a finite time integral and thus the 
integrated autocorrelation time is finite. 

In the critical limit the autocorrelation times diverge. If tr = {T — Tc)/Tc and T^ is the 
critical temperature, for \tr\ ^ we have 

Tint (A;) ~ T^xp,x{f^) ~ ri„t,a;(fc) ~ Itrl"^'' ~ C, (47) 

where u is the usual static exponent and z is a dynamic exponent that depends on the 
considered dynamics: u = 0.683(2) and z = 2.35(2) in the present case.^'^'^^ In the same 
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limit, G{k,t)/G{k, 0) becomes a universal function of the scaling variables 

Q' = ee, S = t/n^,{0), (48) 

i.e. we can write _ 

S^ = r(Q2,S), (49) 

G{k,0) 

where T{Q^, S) is universal, even in S, i.e., T{Q^, S) = T{Q^, —S), and satisfies the normal- 
ization conditions 

r(g^o) = i, / T{o,s)ds = i. (50) 

Jo 
The function G{k, 0) is the static structure factor whose critical behavior has been discussed 
in Sec. Ill A. Using Eq. (15) we can write G{k,t) = xg{Q'^)^{Q'^, S). Analogously, we have 

where the scaling function fmt{Q^) is universal and satisfies /int(O) = 1. 

It is important to note that Eq. (2) does not necessarily imply that the scaling function 
T{Q^, S) decays nonexponentially. On the contrary, as argued in Sec. I, the Griffiths tail (2) 
becomes irrelevant in the critical limit. In view of that discussion it is natural to define a 
scaling function 

/„,,W^).-^ltoj-^, (52) 

which we call, rather loosely, the scaling function associated with the exponential autocor- 
relation time. Indeed, if /exp(Q^) is finite, for S* ^ oo we have 

^(g^5) ~ S'^expi-S/U.iQ')), (53) 

where a is some critical exponent. In terms of quantities that are directly accessible numer- 
ically, we can define it as 

/exp(g^) = hm hm ^^^. (54) 

Of course, the two limits cannot be interchanged. 
The dynamic structure factor G{k,uj) is defined as 

/OO /'OO 

dtG{k,t)e''^* = 2 dtG{k,t) cos cut. (55) 

-OO Jo 
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In the scaling limit we introduce a new scaling function a{Q^,w) defined by 

a(Q^^.) = ^^i^ w^con,,iO). (56) 

T-int(0)G'(A;,0) 

The function o'{Q'^, w) is essentially the ratio of the dynamic and static structure factors 
and is directly related to T^Q"^, S): 

POO 

a{Q\w) = 2 dST{Q\S) cos wS. (57) 

Jo 

It is even in w and satisfies the normalization conditions: 

a(0,0)=2, r'^a{Q\w) = l. (58) 

Moreover, we have (t(Q^,0) = 2/int(Q^)- 

For a Gaussian theory the spin-spin correlation function is given by 

V -n(fc2+r)|t| 

GG{k,t)= ^,^^ . (59) 

It follows Tij^tik) = [r2(A;^ + r)]^^, so that 

^(Q^ s) = e-('3^+m^i, fUQ') = /cxp(g^) = Y^. (60) 

Finally, we have 

ti;2 + (i + Q2)2- 



rn2 ^ 2(1 + g2) 

a(Q ,w) = „ ,, , ^3,3 . (61) 



B. Field-theory results 

The dynamic structure factor can be computed in perturbation theory. The explicit 
one-loop calculation is reported in App. B. Two facts should be noted. First, perturbation 
theory predicts an exponential decay for r{Q'^,S) for any Q^. This is consistent with the 
argument presented in the introduction, which predicted the absence of the Griffiths tail in 
the critical scaling functions. Second, one-loop perturbation theory predicts /cxp(Q^) to be 
independent of Q^. We wish now to argue that this result is exact and is related to the 
breaking of translational invariance in disordered systems. Indeed, consider the spin-spin 
correlation function for a given disorder configuration {p}, 

7(xi, X2; ti - ts; {p}) = {a{xi,ti)a{x2, 12)) p, (62) 
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and the corresponding Fourier transform 

l{ku k2; ti - h; {p}) = J^ e'''''''+"'''''i{x,, x^; h - ^2; {p})- (63) 

In pure systems translational invariance implies that ^{ki,k2;ti — t2', {p}) vanishes unless 
ki = —k2. This is not the case in disordered systems, where translational invariance is lost. 
The average of 7 over disorder vanishes for ki 7^ —k2 [it indeed corresponds to G{k^ ti — 12)], 
and thus translational invariance is somewhat recovered. However, this does not mean that 
the critical theory is translationally invariant. For instance, consider 



\i{k,M-M-t2;{pW ■ (64) 

It can be easily verified in perturbation theory that this quantity is not zero for any ki and 
k2. Note that this breaking of translational invariance survives in the infinite- volume limit 
only close to the critical point. In the paramagnetic phase, far from the critical transition, 
self- averaging occurs and thus also the quantity (64) vanishes for ki 7^ —k2 when L -^ 00. 

Let us now show that, if translational invariance (both for the Hamiltonian and the 
transition rates) holds, the decay rate is k dependent: modes corresponding to different 
momenta decouple. Indeed, following Refs. 9,44, let C be the Liouville operator associated 
with the dynamics, and A^ and ipa be the corresponding eigenvalues and eigenvectors. Then, 
we have the spectral representation (t > 0) 

G{k,t) = J2^~''''MkM^)\', (65) 

a 

where the sum runs over all eigenstates with nonvanishing eigenvalue of C Here we have 
introduced the inner product 

{f\9) = J2''-f*c.9a, (66) 

a 

where / and g are functions defined over the configuration space, tTq, is the equilibrium 
distribution, and the sum runs over all configurations a of the system. If the system is 
translationally invariant, C commutes with the generator T of the translations; hence, the 
eigenstates of C are also eigenstates of T. Thus, we have decoupled sectors corresponding 
to different values of the momentum k and therefore we have 

G(fc,t) = 5]e-^"«*|(a(A;)|^„(A;))p, (67) 
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FIG. 4: (Color online) Scaling function r{Q'^,S) as a function of S, as obtained in one-loop 
perturbation theory. 

where the sum runs over the eigenstates of momentum k. Hence, if Xi{k) is the smallest 
eigenvalue in each sector, we have G{k,t) ~ e"''* with k, = Xi{k); hence, the decay rate is 
k dependent. If translational invariance is lost, all eigenf unctions contribute to each single 
value of k. Note, however, that this does not necessarily imply that the decay rate k in 
Eq. (5) is Q independent. Indeed, one should average over the disorder distribution and this 
average could wash out the effect. We expect this to happen in the infinite- volume limit at 
fixed T, for T ^ Tc. The perturbative results show that this is not the case at the critical 
point. Hence, all modes are coupled in the critical limit and k is momentum independent. 
This argument indicates that the Q-independence of n is strictly related to the breaking of 
self-averaging at the critical point and thus, we expect a similar phenomenon to occur for 
the low-temperature critical dynamical structure factor. 

In Fig. 4 we report r(Q^, S) as obtained by using Eqs. (B29) and (B32) and simply setting 
e = 1. The behavior we observe is quite different from what is observed in the Gaussian 
model. In this case, Eq. (60) implies \nT{Q^, S) = —(1 + Q^)!*?]. As a consequence, with a 
logarithmic vertical scale, the data fall on straight lines with increasing slope as Q^ -^ oo. 
Here instead, T{Q^, S) first decreases rapidly and then bends so that the large- S" decay is 
Q^-independent. This behavior is also very different from that observed in the pure Ising 
model, whose dynamical critical behavior is very close to that of the Gaussian model. ^^ 
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If /cxp(<5^) = /cxp is independent of Q^, for S* ^ oo we expect a behavior of the form 

T{Q\S) ^ /(g2)5"exp(-5//e.p), (68) 

where a is a critical exponent. At one loop, the calculations reported in App. B give a = 
for Q^ = 0, a = — 1 for Q"^ ^ 0, and /(Q^) ~ Q"^ for Q —^ oo. In general, we expect f{Q^) 
to vanish with a nontrivial exponent in the large-Q limit and thus we write 

f{Q') ~ Q^S (69) 

with a new exponent (. 

Given G{k,t), one can compute G{x,t), which can be written in the scaling form 

G{x, t) = ^~d+2^nF(Y\ S), Y^ = x^li\ (70) 

Perturbation theory, see App. B, indicates that FiY'^^S) is not analytic for Y^ -^ 0. It 
predicts a behavior of the form 

F{Y\S) = US) + h{S)\Y\^ + ---, (71) 

where A is a new exponent that can be related to the exponent (^ which appears in Eq. (69): 
\ = C^ — l—rj {md dimensions, as we discuss in App. B, A = C,+2 — d — ri). The exponent A is 
positive (hence C niust be larger than 1 + r^), since G{x = 0, t) is always finite. The quantity 
fx{S)\Y\''^ represents a subleading nonanalytic correction to the leading term fo{S). 

Finally, in Fig. 5 we report the one-loop perturbative expression of o'{Q'^,w). Note that 
the width of o"((5^, w) does not decrease with increasing Q^, as it does in the Gaussian model. 
This is a consequence of the large- S* behavior of r{Q'^,S), whose decay is independent of 

C. Simulation details 

In this section we study the critical dynamics of Hamiltonian (9) at p = 0.8 in the 
high-temperature phase, with the purpose of determining the time-dependent spin-spin cor- 
relation function G{k,t) and the related dynamic structure factor G{k,uj). We perform 
simulations on lattices of size L < 128 in the range 0.275 < (3 < 0.284, corresponding to 
4 < ^ < 16. For each disorder sample, we start from a random configuration, run 1000 
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FIG. 5: (Color online) Scaling function a{Q'^,w) as a function of w = u;Tint(0), as obtained in 
one-loop perturbation theory. 

Swendsen-Wang and 1000 Metropolis iterations for thermalization, and then A^it Metropolis 
sweeps [typically, we took A^it varying between 30rint(0) and lOOrint(O)]. The number of 
samples varies between 5000 and 20000. We measure the second-moment correlation length 
^(/3, L) defined in Eq. (31) and the correlation function G{k,t). As we did for the static 
structure factor, we determine G{k,t) as 



^(^'^) = ^E(^'''' + ^''' + ^''')<^(0'0'0'OM^'^'^''^))' 



(72) 



x,y,z 



where the sum runs over the coordinates (x, y, z) of the lattice sites; the time t is expressed 
in units of Metropolis lattice sweeps. 

Given G{k,t), we determine rint(O). More precisely, we determine rint,x(0) with x = 5, 
as defined by the self-consistent equation (46). As discussed above, this is a good auto- 
correlation time for any x; therefore, we use this quantity to obtain a high-temperature 
estimate of z. We have also determined Tint,a;(0) with x = 8. The results for x = 5 and 
X = 8 are consistent within errors, indicating that we can take Tint, 5(0) as an estimate of 
TintlO). We also consider the effective exponents rexp,a;(0) defined by Eqs. (44) and (45) with 
A(t) = G{k = 0,t). The results we quote correspond to n = 2. 

Some results are reported in Table I. Since we are interested in infinite-volume quantities, 
we must be sure that finite-size effects are negligible. A detailed check is performed at 
(3 = 0.281, where we can compare simulation results at different values of L, corresponding 
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TABLE I: MC results for the randomly site-diluted Ising model at p = 0.8. We report the num- 
ber of samples Ng, the number of Metropolis iterations in equilibrium A'^it, the second- moment 
correlation length ^, the zero- momentum integrated autocorrelation time Tint(O), and the effective 
zero-momentum exponential autocorrelation times rcxp,x(0), for x = 1,2. 



p 


L 


Ns 


Ait 


C 


ri„t(0) 


'rcxp,l(0) 


Tcxp,2(0) 


0.275 


32 


20000 


30000 


4.452(4) 


36.66(27) 


37.92(18) 


39.9(5) 


0.278 


32 


20000 


5000 


5.601(4) 


62.25(23) 


65.70(25) 


71.3(6) 




64 


20000 


5000 


5.622(2) 


61.98(21) 


65.41(30) 


69.4(9) 


0.280 


32 


20000 


8000 


6.872(7) 


102.1(6) 


106.8(5) 


118.2(1.2) 


0.281 


32 


20000 


10000 


7.800(9) 


139.4(7) 


146.7(7) 


164.1(1.7) 




64 


20000 


5000 


7.917(4) 


139.7(8) 


148.2(7) 


158(2) 




128 


10000 


5000 


7.924(2) 


138.8(1.0) 


147.1(1.1) 


157(3) 


0.282 


64 


15000 


20000 


9.331(6) 


207.0(1.0) 


220(2) 


238(3) 




128 


5000 


20000 


9.346(5) 


205.5(1.6) 


218.8(1.3) 


232(4) 


0.283 


64 


20000 


20000 


11.551(10) 


342.8(2.1) 


361.0(1.6) 


402(4) 


0.284 


128 


5000 


50000 


15.837(16) 


716(6) 


753(9) 


842(19) 



to C,{(3, L)/L ^ 0.24, 0.12, and 0.06. No scaling corrections are observed in rint(O) within the 
quoted errors, and thus, for each (3, we assume that the estimate of riiit(O) for the largest 
lattices is an infinite- volume result. Also Tcxp,i(0) apparently does not show finite-size effects. 
On the other hand, Texp,2(0) is clearly decreasing as L increases. This indicates that finite- 
size effects on G{k, t) increase with t, a result that we will check explicitly below, considering 
the correlation function. 

D. Dynamic structure factor 



We first use the estimates of the autocorrelation times to obtain an estimate of z. Since 
the model is approximately improved, ^° the scaling corrections proportional to (/3c — (3)^'^ ., 
ijj = 0.29(3) are suppressed. Thus, the leading scaling corrections behave as {Pc—PY'^^, where 
UJ2 = 0.82(8) is the next-to-leading correction-to-scaling exponent. Hence, r(/9) behaves as 

r(/3) ^ c(/3, - /3)-^^(l + b{(3, - (3)^^ + ■■■), 
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FIG. 6: (Color online) Scaling function G(0, t)/G(0,0) as a function of S, for different values of f3 
and L. 



where A2 = 1^002 = 0.56(6). Thus, we fit the data to 

lnr(/3) = -zu\n{f3, - /?) + a + 6(/5, - (3)^\ 



(74) 



setting (3c = 0.2857431 (3). ^^'^^ If we fit rint(O), including only the data satisfying (3 > /3mm, 
we obtain zu = 1.64(3), 1.59(4), 1.62(8) for /?mm = 0.275, 0.278, 0.280. The resuhs are stable 
with /5min and allows us to estimate zu = 1.61(6) that includes all estimates with their error 
bars. If we now use^° u = 0.683(2), we obtain 



z = 2.36(9), 



(75) 



which is in perfect agreement with the estimate z = 2.35(2) obtained at the criti- 
cal point. ^^ As a check we have repeated the analysis by using Tcxp,i(0). We obtain 
zu = 1.59(3), 1.56(5), 1.47(10) for /?mm = 0.275, 0.278, 0.280, which are essentially con- 
sistent with the estimates obtained above. 

Let us now consider G{k,t). Let us first focus on the case k = Q = 0. Numerical results 
are reported in Fig. 6 vs S* = t/rint(0). Scaling and finite-size corrections are small and 
indeed all data fall approximately onto a single curve. Some deviations are only observed 
for S* > 3, indicating that finite-size corrections increase with S. Let us now consider the 
large-S* behavior and let us estimate the universal ratio /cxp(O). The data show a reasonably 
good exponential behavior so that we can assume that we are considering values of t that 
are much before the region in which G{0,t) shows the Griffiths tail. We perform fits of the 
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mm 



FIG. 7: (Color online) Estimates of [6rint(0)]^^, which converges to /exp(O) in the critical limit. 
The coefficient b is obtained by fitting G(0, t)/G(0, 0), as described in the text, see Eq. (76). In 
fit (a) we consider together the data corresponding to {L = 64, (3 = 0.278), {L = 128, (3 = 0.281), 
and (L = 128, (3 = 0.282). In fit (b) we only consider the resuhs obtained for L = 128, /? = 0.281. 



form 



In 



G{0,t) 



a — bt. 



(76) 



^(0,0) 
including each time only data in the range S^in ^ S < Smax ~ 5. The fit parameter b 

provides an estimate of /exp(O): in the critical limit /exp(O) = [&Tint(0)]^^. Since finite- 
size corrections are important, we only consider data with small ^/L. In Fig. 7 we report 
the results corresponding to two sets of data. In fit (a) we consider three sets of results: 
those corresponding to L = 64, /3 = 0.278 and those with L = 128 and p = 0.281,0.282. 
Correspondingly, we have ^{P,L)/L = 0.088,0.061,0.073, respectively. In fit (b) we only 
use the lattice with the smallest value of ^/L available: L = 128 and f3 = 0.281. The results 
of fit (a) become independent of S'mm for Smm ^ 2.5 and give /cxp(O) = 1.21(1). Fit (b) is 
less stable and a plateau is less evident. They hint at a lower value for the ratio, varying 
between 1.20 (at S^am = 3) and 1.18 (at S^i^i = 4), though with a large statistical error. We 
have also analyzed the data corresponding to lattices with larger C,/L, finding larger values 
of [&rint(0)]^^ This indicates that this quantity decreases with decreasing ^{P,L)/L and 
thus the difference obtained between fits (a) and (b) may be a real finite-size effect. For this 
reason our final result corresponds to fit (b). We quote 



/exp(O) = 1.19(3), 
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FIG. 8: (Color online) Scaling function r{Q'^,S) as computed numerically for two values of Q^. 
For each Q, we report results corresponding to three different values of /3 and L. 

where the error has been chosen conservatively, in order to include the result of fit (a) with 
its error. This result is very close to the one-loop FT estimate. Eq. (B32) gives /exp ~ 1.168 
for e = 1. 

Finally, we provide an interpolation of our numerical data. The curves reported in Fig. 6 
are well fitted by a function of the form 



r(0,^) = e-''^ + 5^afc- 



[cSf 



k=l 



+ (c^)%«^' 



(78) 



The constant n has been fixed by using Eq. (77): we take k = 0.840 ~ 1/1.19. All other 
constants have been obtained by a fit of the data for < 5* < 5. We obtain c = 1.69, ai = 
-0.23825, a2 = 0.16430, as = 0.13261, 04 = -0.28028. As a check of this parametrization 
we verify the normalization conditions (50). The first condition is satisfied exactly, the 
second one to very good precision: the integral between and infinity of r(0, S) as given by 
the parametrization (78) is equal to 1.0033. 

Let us now consider G{k, t) for fc 7^ 0. Again, let us first discuss the finite-size and 
scaling corrections. For this purpose, we must compare G(A;, t) for different values of (3 and 
L, but at the same value of Q = k^- Since the momenta accessible on a finite lattice of 
size L are quantized and therefore estimates are obtained only for Q = 2'n-n^/L, n integer, 
for each t we should interpolate the numerical data as we did in Sec. IIIC However, by a 
fortunate accident, such an interpolation is not needed here. Indeed, the lattice with L = 64, 
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FIG. 9: (Color online) Ratio G{k,t)/G{0,t) ^ r{Q'^,S)/T{0,S) obtained at /? = 0.282, L = 128, 
versus S. Results correspond to different values of Q^ = 47r^n^^^/L^. The values n = 2, 4, 8, 12, 20 
correspond to Q^ « 0.84, 3.37, 13.5, 30.3, 84. 

p = 0.281 and that L = 128, l3 = 0.284 have both ^{(3,L)/L = 0.1237(1). Moreover, for 
L = 128, (3 = 0.281, ^/L is exactly 1/2 (within the small statistical errors) of the previous 
value. Thus, results with the same k for the first two systems and those with 2k for the 
third one correspond quite precisely to the same value oi Q. In Fig. 8 we report results 
corresponding to Q = 27r x 0.1237 ^ 0.78 and Q = 207r x 0.1237 ^ 7.8. All results fall again 
onto a single curve for both values of Q. Finite-size and scaling corrections are apparently 
negligible in this range of values of Q and S. This result should be compared with what 
we observed for the static structure factor in Sec. Ill C. There, a good scaling behavior was 
only observed at fixed Q and not at fixed Q. Here instead, scaling corrections at fixed Q 
are quite small; the behavior at fixed Q is actually slightly worse. 

An important prediction of the FT analysis is that T{Q'^, S) decays with the same rate 
for all values of Q. To check this prediction we consider the ratio G{k,t)/G{0,t), which 
converges to T{Q^, S')/r(0, S) in the scaling limit. For 5 ^ oo, this quantity should behave 
as 

^"exp [SiU^iQY' - /cxp(O)-^)] , (79) 

where a is some critical exponent. Field theory predicts /cxp(<5^) = /exp independent of Q, 
so that we expect T{Q^, 5')/r(0, S) to behave as S*" for large S, without exponential factors. 
Thus, if field theory is correct, hi[G{k,t)/G{0,t)] should become constant as t increases, 
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FIG. 10: (Color online) Ratio Q^G{k,t)/G{0,t) ^ Q^T{Q^,S)/r{0,S) obtained at ^ = 0.282, 
L = 128, versus S. Results correspond to different values of Q^ = Att'^ti'^^'^/L'^. The values 
n = 2, 4, 8, 12, 20 correspond to Q^ « 0.84, 3.37, 13.5, 30.3, 84. The figure on the left corresponds 
to C = 1.9, that on the right to (" = 2.3. 

apart from possible slowly varying logarithmic corrections. In Fig. 9 we show this ratio for 
the lattice with j3 = 0.282, L = 128, which has been chosen because of its relatively small 
errors up to 5 ~ 4. The results for (3 = 0.284, L = 128, which are more asymptotic and give 
access to larger values of Q, are more noisy. The plot shows that the MC data are consistent 
with the FT prediction. Note that the constant behavior is observed better for larger values 
of Q. This is in agreement with the FT results shown in Fig. 4 and can be understood 
quahtatively quite easily. Roughly, at one loop r((5^, S) is the sum of two terms. 



ae 



-s^f,^-ii+Q')s 



(80) 



(we neglect here additional powers of Q and S) , so that the ratio we are considering corre- 
sponds to [a + be~'^ ^)/{(^ + b). Thus, the ratio approaches a constant with corrections of 
order e"*^ '^ . For large Q^ they die out fast, and thus a constant behavior is observed for 
small values of S. 

While the decay rate of T{Q^, S) is independent of Q^, the amplitude decreases rapidly 
with Q^. For large Q^ we expect the behavior r{Q'^,S) ~ S'^Q^^ exp{—K,S), where k = 
1/ fcxp, see Eq. (68). We wish now to obtain a rough estimate of the exponent (. For this 
purpose we take the data that appear in Fig. 9 and we multiply them by Q^, trying to 
fix ( in such a way to obtain a good collapse of the data. In Fig. 9 we report the scaled 
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TABLE II: Numerical values of the coefficients appearing in the interpolation formula (81). 





Q2 = 0.842 


Q2 = 5.262 


Q2 = 21.05 


Q2 = 47.36 


oo 


-2.30388 


0.0972492 


-0.1390160 


-0.1482217 


ai 


1.85440 


-0.1867741 


-0.0166780 


0.0609449 


02 


-1.03426 


0.2245075 


0.0882145 


-0.0108184 


as 


-1.04434 


-0.0670895 


-0.0483811 


-0.0040681 


04 


2.29488 


-0.0788428 


0.1494524 


0.1510844 


do 


3.30388 


0.9027518 


1.1390167 


1.1482211 


di 


-1.87228 


-0.5008092 


-13.808512 


-19.439356 


d2 


2.03391 


12.846119 


114.32824 


182.53584 


ds 


-0.64676 


-27.878988 


-319.65479 


-587.00842 


d4 


0.13949 


58.784027 


519.97023 


1021.6385 


K2 


1.95 


9.08 


10.46 


12.40 


C 


2.000 


1.124 


1.087 


1.282 



results corresponding to two different values of (. If we try to have a good collapse of the 
data corresponding to n = 4, 8, 12 the best result is obtained for ( = 2.3. However, the 
data with n = 20 behave in a significantly different way. If we try to include also the data 
with n = 20, the quality of the collapse worsens and the best result is obtained for ( = 1.9. 
These results indicate that C ~ 2 (but with a large error), so that T{Q^, S) behaves roughly 
as Q^'^S"' exp{—K,S). It is interesting to observe that this is exactly the behavior predicted 
close to four dimensions by perturbation theory. 

Finally, we determine an interpolation formula for r{Q'^,S). We find that all data are 
well fitted by taking 

- , (cSY- 

k=l ^ ' k=0 

fixing K,i = 0.840. The results of the fits for a few chosen values of Q^ are reported in Table 



^(Q^ S) = aoe-^ + J^ ^^ i,LL. + E ^^^'^""'' 



(81) 



II. We have not required a^ + do = 1, a condition that follows from r((5^, S* = 0) = 1, but 
we have verified that the results satisfy this condition quite precisely. By using a linear 
interpolation, the results we report should allow the reader to determine T{Q'^, S) for any 
Q in the range < Q^ < 50 with reasonable precision. We stress that this interpolation 
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FIG. 11: (Color online) Numerical estimate of the scaling function a{Q'^,w) as a function of 
w = a;Tint(0), as obtained by integrating the interpolating expression (81). The values of the 
momenta are: Qo = 0, Qf = 0.842, Q^ = 5.262, Q| = 21.05, Ql = 47.36. 

formula only represents a compact expression that summarizes the numerical results. The 
chosen parametrization has indeed a purely phenomenological value. 

From these expressions it is easy to determine the scaling function (j{Q^, w) related to the 
dynamic structure factor. In Fig. 11 we plot the scaling function a{Q'^,w) as obtained by 
integrating the interpolating function determined above. We report it for the same values 
of Q^ that appear in Table II. The qualitative behavior is in full agreement with the FT 
prediction, compare with Fig. 5. Quantitatively, perturbation theory is also reasonably 
predictive. For Q = and w < 5, the relative differences between the FT and the numerical 
expression are less than 2%. For larger values of Q^ differences increase: for instance, for the 
values of Q that appear in Fig. 11, field theory predicts cr{Q^, 0) ^ 2, 1.09, 0.32, 0.091, 0.041, 
while we obtain numerically a{Q'^,0) ~ 2, 1.09, 0.33, 0.105, 0.072. These discrepancies are 
probably the fault of both field theory — after all, we are at one loop — and of the numerical 
results — for large Q the data are noisy and the estimates have a large error. In any case 
this comparison indicates that, up to Q = 5, errors are under control and the reported 
expressions are precise enough for all practical purposes. 
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APPENDIX A: PERTURBATIVE RESULTS FOR THE STATIC STRUCTURE 
FACTOR 

The static behavior of Ising systems with random dilution can be studied starting from 
the Landau- Ginzburg- Wilson Hamiltonian*^^ 



n= (fx 



■^{d^^^f + -rv?^ + -V^v^^ + ^^oV^^ 



(Al) 



^(^) = /- exp 

yAixw 



(A2) 



where ^(x) is a spatially uncorrelated random field with Gaussian distribution 

Aw_ 

Using the standard replica trick, it is possible to replace the quenched average with an 
annealed one. As a result of this procedure, one can investigate the static critical behavior 
of RDI systems by applying standard FT methods to the Hamiltonian*^^ 

Hreplica = / d''x{j2 \ [(^M^^)' + ^0'] + E ^i ("^ + ^0%) 0-0?}, (A3) 

i ij 

where i,j = 1, ...N and uq = —6w. RDI results are obtained by taking the limit A^ ^ 0. 

1. y^-expansion results 

The scaling function g{Q'^) can be determined by using the results reported in Ref. 37. 
We obtain the expansion 

9{y)-' = l + y-4^2(y) + 18^^[24 + 7C(3)]63/^^2(y)+0(62) 

= l + y- 0.0188679e^2(l/) + 0.0698876^/^^2(1/) + ^(e^) , (A4) 

where ip2{y) is the two-loop contribution defined in Ref. 37. Note that the only relevant three- 
loop diagram contributes only at order e^/^ (hence, at five loops), since it is proportional 
to 

^{u + 3v){4u + 3vf, (A5) 

and, at the fixed point, 4u* + 3v* is of order e and not of order y^. 

The expansion of ip2{y) for small y can be found in Ref. 37. It allows us to obtain the 
expansions 

C2 = 0.000141891e- 0.0005255676^/2 + 0(e2), 
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C3 = -3.62134 ■10"^e + 0.0000134135e^/2^O(e2), 

C4 = 1.53623 -lO'^e- 5.69021 ■10"V/^ + O(e2), 

C5 = -8.28575 ■10"^e + 3.06905 ■10^V/2 + O(e2). (A6) 



The expansion of ip2{y) for large values of y is^^ 

My) = -^y logy + 2goy - ^ log' y + 2Q, + ... (A7) 

with Qo ^ 0.507826 and Qi ^ 0.1289. Matching the large-momentum expansion of Eq. (A4) 
with the Fisher-Langer behavior (20) we obtain the expansions of the coefficients Cf 

Ci = l + 0.0191632e - 0.07098096^/^ ^ q(^^2^^ 

C2 = -1/2 - 0.757042e^/' + 1.34297e + 02,36^^' + 0{e^), 

C3 = -1/2 + 0.7570426^/^ - 1.35726e + 03,36^/=^ + ©(e^), (A8) 

where 02,3 + 03,3 = 0.052964. Setting e = 1, we obtain Ci ^ 0.95, ^2 + ^3 = -0.96. 

2. Massive zero-momentum results 

We have determined the low-momentum behavior oi g{Q'^)~^ in the MZM scheme by using 
the perturbative results of Ref. 37. The four-loop expansions of the first few coefficients c„ 
are: 

C2 = —u^ —uv —v^ - 0.0000120404^3 - 0.0000481617^2^ 

6480 2430 10935 

-0.0000481617Mt;2 - 0.0000142701^;^ - 0.00000718972m^ - 0.0000383452^3^ 

-0.0000617074mV - 0.0000379116Mt;3 - 0.00000842479^;^ 

C3 = — ^ — u^ + -^^uv + v^ - 0.00000281924^3 - 0.00001 12769^2^; 

122472 45927 413343 

-0.0000112769Mt;=^ - 0.00000334132t;3 + 0.00000660026m^ + 0.00000352014^3^; 

+0.0000056289lM2t;2 + 0.0000033937lMt;3 + 0.000000754158^;^ 

111 

C4 = u"^ - UV v^ + 0.000000347995^3 + 0.00000139198^'^; 

1889568 708588 1594323 

+0.00000139198Mt;2 + 0.000000412438t;3 - 0.000000141114^^ - 0.000000752609M3t; 
-0.00000119627M't;' - 0.000000740733Mt;3 - 0.000000164607t;^ (A9) 

The MZM renormalized quartic couplings u and v are normalized so that at tree level 
u = Uq/itl and v = Vq/itl. Their fixed-point values are u* = —18.6(3) and v* = 43.3(2) 
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(obtained by means of MC simulations^'^), and u* = —13.5(1.8) and v* = 38.0(1.5) (obtained 
by resumming the six-loop /3-function^^). 

APPENDIX B: ONE-LOOP CALCULATION OF THE RESPONSE AND COR- 
RELATION FUNCTIONS 

The relaxational model-A dynamics is described by the stochastic Langevin equation^ 

where ip{r,t) is the order parameter, 7i((/9) is the Hamiltonian (Al), ^2 is a transport coeffi- 
cient, and C{r,t) is a Gaussian random field (white noise) with correlations 

(C(r, t)) = 0, (C(ri, ti)C(r2, h)) = n6{n - r2)6{t, - t,). (B2) 

The correlation functions generated by the Langevin equation (Bl) at equilibrium, averaged 
over the noise ( and the quenched disorder ip, can be obtained from the FT action^^ 

S{ip, 0) = / d'^x / dt ifidtif - nAip -nip + nrip) (B3) 

where ip is the response field. In this framework, no replicas are introduced.'*^ We consider 
the correlation function G{x,t) and the response function R{x,t) defined by 

G{X2 -xi,t2 -ti) = {ip{Xi,ti)(p{x2,t2)), (B4) 

R{x2 -xi,t2- h) = {ip{xi,ti)ip{x2, t2)), (B5) 

and their spatial Fourier transforms G{k,t) and R{k,t). In equilibrium, they are not in- 
dependent, but related by the fiuctuation-dissipation theorem; for t > they satisfy the 
relation dtG{x,t) = —QR{x,t). For a general introduction to the FT approach to equilib- 
rium critical dynamics, see, e.g., Refs. 47,48. Some perturbative calculations can be also 
found in Refs. 49. 



1. One-loop calculation 

We first compute the response function R{k, t) and then use the fiuctuation-dissipation 
theorem to derive the correlation function G{k,t). At one-loop we obtain in dimensional 
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FIG. 12: The one-loop graphs contributing to the response function. We indicate with Lg{k,t) 
and with Lu{k,t) the contribution of the graph on the left and on the right, respectively. 

regularization 

R{k, t) = Rcik, t) - ^Lg{k, t) - y L„(A;, t) + 0{ul u^go, gl), (B6) 

where 

Rcik, t) = e{t) exp[-fi(P + r)t\ (B7) 

is the Gaussian tree-level response function, and Lu{k,t) and Lg{k,t) are the one-loop con- 
tributions, see Fig. 12. It is straightforward to obtain 

Lg{k, t) = nl dsl —^ Rcik, t - s)RGik, s)GGiq, 0) 

= _!^r^--/^ntRG{k,t), (B8) 

where e = A — d, Nii = 2/[F((i/2)(47r)'^/^], and Gcik^t) is the Gaussian tree-level correlation 
function (59). Analogously, we obtain 

Lu{k, t) = n^ dsi / ds2 / — — ^ Rcik, t - 52)^0(9, -52 - Si)RG{k, Si) 

= —{l + e^E/2){nty/\nte-l)RG{k,t) + e{t)-^e-^''F{nte)+0{e), (B9) 



(47r, 



where 

F{x) = -1 + e-^ + e'^'ix - l)[Ei(a;) - 7^ - Inx], (BIO) 

and Ei(x) is the exponential integral function. 

Renormalizing the response function in the MS scheme we obtain 

R{k, t) = RG{k, t) - ^^tr \nrRG{k, t) (Bll) 



-{-lE + \ri^t){^te - l)RGik,t) - -—9it)e-''''Fintk'), (B12) 
487r^ 487r^ 

where r, Q, u, g are renormalized parameters (note that we use the same symbols r and Q 
for both the bare and the renormalized parameters, since no confusion can arise). The final 
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expression is obtained by setting g and u equal to their fixed-point values: 



/ = 327r^/^, 



u 



., ;6e 

-247r^/— . 
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The correlation function G{k, t) for t > can be obtained from 

G{k,t) = n R{k,t), 



(B13) 



(B14) 



which follows from the fluctuation-dissipation theorem and the fact that G{k, t) ^ as 
t ^ oo. A straightforward calculation gives (in the following we always assume t > 0) 



G{k,t) 
G(A;,0) 



-nt(fc2+r) 



g 



u 



487r2' 


W ^ i 


u 

487r2 


ntk'^ 
\ 


u 


r 



327r2 

r 



2 -Qt(fc2+r) _ 



U 



-nt{k'^+r) _ -Q.tr 



k"^ + r 



487r2 



where 



487r2 A;2 + r 



^(^,0) 



Ei(-l]tr) - e 



-nt(fc2+r) 



(7£; + In fitr) 



1 



1 , , , rlnr 

-{g + 2u/?,)- 



(B15) 



(B16) 



A;2 + r 327r2'^ ' ^(p + r)2' 

It is easy to derive the critical correlation function at one-loop order from Eq. (B15). Taking 
the limit r ^ 0, we obtain 



G{k,t) 
G(A;,0) 



^.ntk^ _ ^(^^ + In nmtk^-'"'^' "" 



487r2 



487r2 



-ntk'^ _ 1 



487r^ 



rttk^e'^"" {Ei{Qtk^) --fE- In QtP). 



(BIT) 



In the critical limit the correlation function should only depend on the scaling variable 

K^ = A;2(r2t)2/^. Keeping into account that z = 2 — u/{2A'k'^) at one loop, we obtain at this 
order (we set u = u*) 



G{k,t) 



e-^ - 



u^ 



487r2 



-K^ 



1 - 



M' 



487r2 



K^e'^' (Ei(i^') - In K^) . (B18) 



G(A;,0) 

This scaling function has a regular expansion in powers of K^ for i^ — > 0, while, for K —^ oo, 
it behaves as 



G{k,t) 
G(A;,0) 



u^ 



1 



'6e 



487r2ir2 2K'^\ h?,' 



(B19) 
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Given Eq. (B18), we can compute G{x,t) at the critical point. We expect the scahng 
behavior 

G{x,t) = {nt)-^'^+"~^^/'F{X), X = x{nt)-^/\ (B20) 

Using Eq. (B18) we obtain 

nX) = 4^(1 - e-^^/^) - ^.FnooAn (B21) 

with 

1 roc 

i^iioop(X) = ^^^ J dK Ji{KX) [e-^' - 1 + K'e-''" {Ei{K^) - IniT^)] , (B22) 

where Ji (x) is a Bessel function. In the derivation we have taken into account that u ~ ^/e. 
It is interesting to note that -Fiioop(-^) is not regular for X ^ 0. Indeed, the explicit 
calculation gives 

i^iioopW = -^ln^ + 0(X2lnX2). (B23) 

This result indicates that F{X) is not regular for X ^ 0, but behaves as 

F(X)^ao + 6o|^|^ + --- (B24) 

where A is a critical exponent. Comparing this expression with Eqs. (B21) and (B23) we 
obtain 

ao + bo = -^ + OiV~e), (B25) 

It is not possible to compute A, Qq, and Bq separately at this order. A two-loop computation 
of the term proportional to e In X is needed. 

In the high-temperature phase, it is convenient to replace r and ^2 with the correlation 
length ^ and the zero-momentum integrated autocorrelation time Ti^t{k = 0), defined in 
Sec. Ill A and IV A, respectively. A tedious calculation gives 

r' = r + ^{g + 2u/3)rlnr, (B27) 

_ 1 g rlnr u fc^lnr u 1 mos^i 

^inH ) - -^^T^ - 327^2 (F + rf ^ 487r2 (P + rf ^ AStt^ k^ + r' ^ ^ 
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Using these results, we obtain the one-loop expression of the scaling function T{Q^, S): 



r(g2,5) 



u 



5(Q^+i) _ " 5(g2 ^ i)e-5(Q^+i) 



487r2 



u 



AStt' 
u* 



487r2 

u* 1 
"487r2Q2_ 



,(7i? + ln5)5g'e 



Q'S- 



2^-SiQ'^+l) 



U 



g-5(Q2+i) _ ^-s 



487r2 

-5(Q^ + l)(Ei(5g2)_^^_l^^g2) 



1 L 



g2 + i 



(B29) 



where -u* is the fixed-point value (B13) of u. 

It is interesting to discuss the large-S* and small- S* behavior of the scaling function (B29). 
For large S we obtain 

u* i + g^e-^ 



T{Q\S)^e 



-Sil+Q"-) 



W 



'l + Q'-Q'\nQ')Se 



~S{1+Q') 



(B30) 



487r2 (g2)2 52 48^2' 

For g^ = the dominant term is the last one. In this case we can rewrite the large-S* 
behavior as 



r(g2, s) ^ exp 



-SI 



u 



487r2 



which gives for the exponential autocorrelation-time scaling function [see Eq. (52)] 



/exp(O) = 1 - 



U 



487r2 



l + -\- + 0{e). 
2 V 53 ^ ' 



(B31) 



(B32) 



For g2 ^ 0, the dominant term is the second one, so that for any Q^ the scaling function 



decays as e / S . Thus, at this perturbative order, we obtain the result 

/exp(g') = l+0(v^). 



(B33) 



The correlation function T{Q'^^S) decays with the same rate for all values of Q. As we 
discuss in Sec. IV B, this is a consequence of the loss of translational invariance in dilute 
systems. 

Equation (B33) should be contrasted with the result obtained for the integrated autocor- 
relation times. Using Eq. (B28), we obtain 



Tint (A;) 



fUQ' 



(B34) 



rint(O) ^"'^^" ^ g2 + l' 

without one-loop corrections. This shows that Tint (A;) decreases as Q increases as it does in 
the Gaussian model. 
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Let us now consider the limit S* ^ 0. The scahng function has an expansion of the form 

T{Q^,S) = l + J2S'"{an + bn\nS). (B35) 

n=l 

Note the presence of terms proportional to In S. They should be generically expected, since 
in the critical limit (it corresponds to Q ^ oo and S ^ 0), the correlation function depends 
on 

The presence of these logarithms implies that the function r((5^, 5*) is not analytic for S* ^ 0. 
It is also important to discuss the large-momentum behavior of T{Q'^, S). For Q'^ ^ oo 
the tree-level term vanishes exponentially as e"'^'^ , while the one-loop term decays only 
algebraically, as l/Q'^- More precisely, for Q^ ^ oo we have 

^(Q^ S) ^ -^ (^^ + Ei(-5)) Q-\ (B37) 

The presence of these slowly decaying terms implies the singularity of the behavior of G{x, t) 
for X ^ and for any t. In the critical limit we expect the scaling behavior (70), i.e., 

G{x,t) = ^'d+2-rip(j2^ g^^ ^j^l^ y2 ^ ^2/^2^ ^g ^^^g^j^^ for F ^ 

F{Y^ ^) - Y^ (^ + Ei(-5)) (l + ^ In >^) + ^*fcoAS) + ... (B38) 

for Y ^ 0, where fcorr{S) is a function of S. The presence of a term proportional to 
InY implies that F(Y'^,S) is not analytic as F ^ 0, i.e. has a behavior of the form 
F(Y^, S) = /o(S') + /a(S')|F|'^ + - ■ ■, where A is the same exponent that appears in Eq. (B24). 
Note that Eq. (B38) is apparently consistent with the assumption that /o(>S') = 0. If this 
were the case, we would obtain 

\ = ^ = -\ — + 0(e). (B39) 

247r2 V 53 ^ ^ ^ ' 

This result would imply A < and thus F{Y'^, S) would diverge as F ^ for any S, at 
least for e small. This behavior is clearly unphysical; thus, /o(>S') should be nonvanishing. 

The critical limit is obtained by taking S ^ 0. Requiring the limiting function to be of 
the form (B20), we obtain 

fx{S) ~ 5(-A+2-^~d)/^^ ^34Q) 
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for S* ^ 0. The S'-dependent prefactor appearing in Eq. (B38) behaves as 1/ S for S* ^ 0, 
which is consistent with these expressions. 

The nonanalytic behavior of FiY"^^ S) a,sY ^ 0, imphes that T{Q'^, S) should decay as 
a power of Q as Q ^ oo. A simple calculation gives 

V{Q\S)^f\{S)Q''-''~^-\ (B41) 

for Q^ — ^ oo. The exponent C, defined in Eq. (69) is given by 

( = A + rf + r7-2 = 2 + 0(v^). (B42) 

Finally, we report (j{Q'^^w)^ cf. Eq. (56). A long calculation gives 
,^n , 2a u* \ aia^ — w'^ + 2aw'^) , aia^ — w"^ — 2a) ^ ,_ n. 



a^ + w^ 247r2 [ w{a^ + w^f 2{a'' + w^f 

w-' + l {w^ + l){a'' + w'') {a^ + w^fy ^ ' 

where a = 1 + Q^. 

For large w, a{Q'^,w) behaves as 

a(g^w)^ — 



u 



(B44) 



which is compatible with the expected behavior ■u;~(2~''+^)/^. 

Note also that the singularities of a{Q'^^w) in the complex w-plane that are closest to 
the origin are w = ±i, independently of Q^. This is a direct consequence of the fact we have 
already noticed that the large-t behavior is momentum independent. As a consequence, the 
width of the structure factor does not decrease with Q^ as it does in pure systems. 
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